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1.  Introduction 

It  is  known  (e.g.  see  Brillinger  (1974))  that  certain 
suitably  defined  estimates  of  the  spectral  density  matrix 
of  the  Gaussian,  stationary  multivariate  time  series  are 
approximately  distributed  as  complex  Wishart  matrices.  So, 
the  problems  of  inference  on  the  covariance  matrices  of  the 
complex  multivariate  normal  distributions  are  closely  re¬ 
lated  to  the  problems  of  inference  on  the  spectral  density 
matrices.  For  a  review  of  the  literature  on  complex  multi¬ 
variate  distributions  and  their  applications  in  time  series, 
the  reader  is  referred  to  Krishnaiah  (1976). 

Motivated  by  the  applications  in  the  area  of  inference 
on  multiple  time  series,  we  investigate  asymptotic  expressions 
in  the  null  and  nonnull  cases  for  the  distribution  of  certain 
power  of  the  likelihood  ratio  statistic  for  testing  the 
hypothesis  that  the  variables  are  independent  and  have  a  com¬ 
mon  variance.  These  expressions  are  in  terms  of  beta  series. 

In  the  case  of  null  distribution,  it  is  found  that  the  ac¬ 
curacy  of  the  approximation  by  taking  the  first  term  alone 
in  the  asymptotic  series  is  sufficient  for  practical  purposes. 
Here,  we  note  that  Krishnaiah,  Lee  and  Chang  (1976)  approxi¬ 
mated  the  null  distribution  of  certain  power  of  the  likeli¬ 
hood  ratio  test  statistic  for  sphericity  with  Pearson's  Type 
I  distribution.  But  this  approximation  is  based  upon  empirical 
study.  In  the  analogous  real  case,  Khatri  and  Srivastava  (1974) 
derived  the  nonnull  asymptotic  distribution  in  terms  of  chi- 
square  series.  In  the  final  section  of  this  paper,  we  discuss 


raas2 


smsasa&ai 
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the  applications  of  our  results  to  the  area  of  inference  on 
multiple  time  series. 
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2.  Preliminaries 

In  this  section,  we  define  some  notation  and  give  some 
lemmas  which  are  needed  in  the  sequel. 

I 

The  Mellin’s  integral  transform  of  a  function  f(x)  of 
real  variable  x  defined  for  x  >  0  is 


f  °° 

M{  f (  •  )  I  t  }  «  XC  1  f (x) dx 

j  0 


(2.1) 


where  t  is  a  complex  variate  (Titchmarch  (1937)). 

c  —  3. 

Lemma  1 .  Let  y  f(y)  be  a  measurable  function  in  (0,®) 
and 

r“ 

f  (y )  dy  <  “>  . 


r°° 

f  c-l 
y 

J0 


Also  let  f  (y)  be  of  bounded  variation  in  the  neighborhood 
of  the  point  y  *  x.  Then 


c+i® 

r 


hi  f (x+0)  +  f  (x-0)  }  = 


2iri 


M{  f  (  •  )  |  t  }  x  t  dt 


c-i® 


where  M{  f  (  ■  )  |  t  }  for  t  =  c  +  iv  exists. 

In  the  .sequel,  we  shall  assume  that  f(x)  is  absolutely 
continuous  in  the  interval  (0,1).  Hence  f(x)  is  of  bounded 
variation  in  the  neighborhood  of  x  of  interest.  Furthermore, 

c+i“ 

f(x)  =  —j-  |  M{  f  ( • )  j  t }  x_t  dt.  (2.2) 


c-  ic 

Note  that  when  f(x)  =  (1-x) 

rl 


b-1 


0  <  x  <  1 .  Then 


i  t-1  ,,  sb-l  ,  ”(t)T(b) 

M{  f  (  * )  I  t }  =  [  x  (1-x)  dx  *=  -  r 


(2.3) 
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r  (a)  *  7rP(p-l)/2  P  r  (a-i+1)  . 
P  i-1 
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The  complex  multivariate  hypergeometric  coefficient  is 


given  by 


where 


H  (  a  -i  +  1 ) 


(a)k  =  a(a+l) . . . (a+k-1) 

k  =  (k.  , . . .  ,k  )  is  a  partition  of  the  integer  k  such  that 
1  P 

k..  >...>  k  >0  and  k  =  k,  +  ...+  k  .  The  transpose  and 

l  —  —  p  —  1  p 

conjugate  of  a  complex  matrix  B  are  denoted  by  B1  and  B 

respectively.  Also,  let  C^(A)  denote  the  zonal  polynomial 

» 

of  a  Hermitian  matrix  A,  (i.e.,  A  =  .  ).  In  addition, 

~  “  _  [ 3 1  ^  c  *  *  •  f  a*.  ^  k-  C  (A) 

rFq(a1»  •  •  •  *ar;  bL . b  5  A)  =  l  l  fb-_ )  ..  .  [b  ]  kT 

k=0  K  I  k  (J  K 

^F^(a^,.,.,a^j  b^,,,.,b^j  A ,  B) 

«  [a.]  ...[a  ]  C  (A)  C  (B) 

y  y  l  ie  r  c  k  _ < 

k-0  k  ^l^K***^q^K  C  <1  )  k' 

M  K  p 

where  a^,...,a^,  b^,...,b^  are  real  or  complex  constants. 
Throughout  this  paper,  etr  B  denotes  the  exponential  of  the 
trace  of  B. 

Lemma  4.  Let  £ :  pxp  be  a  Hermitian  matrix. 


Then 


0V-E'1,  L>  -  «'trL/q  „F0(i  M.L) 


(2.8) 


where  M  =  I  -  qZ  ,  L:  pxp  is  a  positive  definite 
Hermitian  matrix,  q  is  a  constant  and 


0F0(B’D)  =  l  l  JS. 


ck(b)  C^CD) 


k=o  k  c  (i  )  k: 

<  P 


(2.9) 


Proof :  By  using  the  splitting  formula  (James  (1964)),  we 

know  that 

ULU ' , I ) dU  (2.10) 


u(p) 


where 


UlBM)  -  e*" 

dU  is  the  invariant  Haar  measure  on  the  unitary  group  U(p) 
normalized  to  make  the  total  measure  unity.  Since 

.tr  I_1ULU’  =  —  tr  L  -  tr(-  MULU  ’  )  , 


.-1 


(2.11) 

where  M  =  I  -  q£  J'.  We  prove  the  result  by  using  Eq .  (2.11) 
in  Eq.  (2.10) . 

Lemma  5 .  Let  V  be  a  pxp  positive  definite  Hermitian  matrix, 
and  T:  p*p  be  an  arbitrary  complex  symmetric  matrix.  Then 


exp(-trV/q) |v|b_P  (tr  V)m  (TV)dV 
V'  =V>0 

=  f  ( b  ,  < )  r  (bp+k+m)  C  (T)  qbp+k+tn/r  (bp+k)  (2.12) 

P  * 

where 

f  (b ,  k  )  =  ^pCp-1)  n  r(b-i+l+k.),  k  =  (k,,...,k  ). 

P  i=l  1  1  P 

Proof .  We  know  (see  Khatri  (1  966))  that 

f  exp (-tr  V )  | V | b  "P  C^(VT)dV  =  Tp(b,K)  (T) . 
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l<  Jr 

Substituting  Z^VZ  "  for  V  in  the  above  equation  with  Jacobian 
|z|P  (Khatri  (1965))  where  Z^  is  a  Hermitian  positive 


definite  matrix,  we  have 


.  1  ■ 
V' =V>0 


exp (-tr  ZV )  1 V I  P  C  (VT)dV  =  r  (b.x)  C  (TZ  )!z 

K  p  K 


Now  take  Z  =  (—  -  x) 

9 

f  exp (-tr (i  -  x)V)|v|b_p  C  (VT)dV  =  T  (b,ie)  C  (T)(i  -  x)-<bP+k) 
^  J  4  P  K  q 

V' -V>0 


Equating  the  coefficient  of  x  /ml  for  both  sides  of  the  above 
equation  we  obtain  (2.12). 

Lemma  6.  For  any  integer  r,  variate  x  and  Hermitian  positive 
definite  V,  we  have 


“  x  C  (V) 

£  £  - - -  =  x  (tr  V) r  etr  (xV) 

k=r  k  (k-r ) l 


(2.13) 


00  x  a.  (  k  )  C  (V) 

l  i  — i 

k  =  0  k  k! 


=  (x2tr  V2-x  tr  V)  etr(xV)  (2.14) 


03  x  a .  ( < )  C  (V)  .. 

V  r  1  <  r  r  +  2^  TT2,.  ITNr  r+1  , .  „vr+l 

L  l  -  =  {x  tr  V  (tr  V)  -x  (tr  V) 

k=r  k  (k-r)I 

+  2r  xr+1  tr  V2(tr  V)r  1-rxr(tr  V)r 


+  r(r-l)xrtr  V2(tr  V)r  2}etr(xV) 


(2.15) 


«  x  a.  (<)C  (V)  ,  _  „  .  to  o 

£  £  — - - -  =  {x  (tr  V  )  +  4x  tr  V  -2x2  tr  V  tr  V2 

k=0  k  k! 
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+  3x2(tr  V)“-4x2  tr  V2  +  x  tr  V}  etr(xV) 


“  xka  (<)C  (V)  _  -  9  0  ~ 

I  I  - - - - -  =  {  2x 3  tr  V3  +  3x2(tr  V)2-3x2tr  V2 


(2.16) 


k  =  0  <  k ! 


where 


+  2x  tr  V}  etr (xV) 


r 

‘ ,  (  *  )  =  I  k  ( k  -  2j  ) 
j  =  l  2  2 


a~  (<)  =  2  l  k  (k2-3jk  +3j  2) 
3=1  3  3  3 


(2.17) 


(2.18) 


and  k  was  defined  earlier, 


The  above  lemma  was  proved  by  Hayakawa  (1972). 
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3-  Asymptotic  Null  Distribution  of  the 
Likelihood  Ratio  Test  Statistic  f oi 
Spher icity 

In  this  section,  we  derive  an  asymptotic  expression 
for  the  null  distribution  of  the  likelihood  ratio  statistic 
for  testing  the  hypothesis  of  sphericity  for  complex 
multivariate  normal  distribution.  The  expression  obtained 
is  in  the  form  of  a  beta  series. 

Let  Z:  pxl  be  distributed  as  a  complex  multivariate 
normal  with  mean  vector  u  and  covariance  matrix  Z.  The 
density  of  Z  in  this  case  is  known  (see  Wooding  (1956)) 
to  be 

£(Z)  -  A -  stt [-Z-1 (Z-u) cz-u)  '  ] .  (3.1) 

"PUI  -  -  -  - 

Next,  let  Z^,...,Z^  be  N  independent  observations  on  Z 
and  let 

N  _  - 

A  -  l  (Z  -Z.) (Z  -Z.)  =  (A  )  (3.2) 

t  =  l  ~ c  ~ 

where  Z  denotes  the  conjugate  of  Z,  and 

N 

NZ.=  l  Z 
t  =  l 

We  are  interested  in  testing  the  hypothesis  H  where 
2  2 

H:  Z  =  a  I  and  a  is  unknown.  The  hypothesis  H  can  be 
P 

decomposed  as  H  =  A  where  is  the  hypothesis  that  Z 
is  diagonal  matrix  and  H2  is  the  hypothesis  that  the  diagonal 
elements  of  Z  are  equal  given  is  true.  The  likeli¬ 
hood  ratio  test  statistic  for  testing  and  H0  are  known  to 


be  A^  and  A2  respectively  where 


pr- 


X2  = 


P  ^ 

11  Aii 
i=l 


P 

11  A?i 
i-1 

(trA/p)pn 


(3.3) 


(3.4) 


where  n  =  N-l .  The  likelihood  ratio  test  for  H  is 


1  / 

Now,  let  w  =  \  n.  Then 


x  =  x1x2  . 


P 

n  i 


a  v 


(3.5) 


(3.6) 


where  _>.  .  .  _>  l  are  the  eigenvalues  of  A.  The  moments 


of  w  under  the  hypothesis  H  are  known  to  be 

E  (wh)  =  [pph /Tp (n) ] T (np)Tp (n+h ) /T (np+ph) 


(3.7) 


Next,  let  u  =  w  where  w  is  given  by  Eq .  (3.6)  and  s  is  a 

constant  to  be  chosen  to  govern  the  rate  of  convergence  for  the 

t  h 

resultant  series.  The  null  h  moment  of  u  is  obtained  by 
replacing  h  with  h/s  in  Eq .  (3. 7). By  using  the  Mellin's  in¬ 

version  transform  (see  Eq.  (2.2)),  the  density  of  u  becomes 


c  +  i®  n  T(n  +  -  i  +  1) 

.  f  u-»-i  pph/s  ±=l - L_ - dh 

2ui  cJ_iM  T(np  +  P  ~) 


(3.8) 
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and  K(p,n)  =  T(np)/  [  II  T(n  -  i  +  1)].  Set  m  =  n  -  5, 

i  =  l 
Vi  t 

d  =  c  +  ms  and  m  +  —  =  —  ,  where  6  is  also  a  converging 
factor  to  be  chosen  for  the  resultant  series.  Then,  we  have 


.  /  \  ^  /  \  ™  p  tn  s  m  “1  1 

f(u)  -  K ( p , n )  p  u  - 


u  C  <j)  (  t)  dt 


(3.9) 


d-ico 


n  r(-+6-i+i) 


*(t)  »  ppt/s 


(3.10) 


r(pf  +  6P) 


By  the  use  of  the  formula  for  the  asymptotic  expansion  of 
gamma  function 

log  T  (x  +  b)  =  log  /2tt  +  (x  +  b  -  y)  log  x  -  x 


-  I  (-D 


t  Br+l(b> 

r  (r  +  1)  x1 


(3.11) 


for  b  bounded  and  Br(b)  is  the  Bernoulli  polynomial  of  degree 
r.  So,  we  have 

Ezi  :  a,  A 

log  <j>(t)  =  1  og  ( 2  it  )  *  sv  p-<sP+i5  +  log  t  v  +  l  -£•  (3.12) 


r  =  l  t 


where 


v  =  (p  -l)/2 


(3.13) 


(  , nT  r  B  .  ( 6  p )  P 

Ar-  r(hlt  ^ - J,  <3’14) 


Hence 


t(.)  -  (2,)  2  sv  p-6p+li 
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x  t 


-v 


[1  +  I  -f  3 

r  =  l  tr 


(3.15) 


The  coefficient  can  be  obtained  by  the  recursive  equation 
Eq.  (2.7), 


«r  *  7  tAt  Vt  •  %  ’  »• 


(3.16) 


Since  ( t )  =  0(t  V)  ,  we  can  write  $(t)  as  follows  by  applying 
Lemma  2 : 


-v 


00  Q  “  T(t+a) 

fi  +  i  -f>  -  i  *4 - 

r =1  t  i=0  T (t+a+v+i) 


(3.17) 


and  a  is  a  constant  to  be  chosen  to  govern  the  rate  of  con¬ 
vergence  for  the  resultant  series.  Using  Eq .  (3.11)  to  expand 

the  gamma  function  on  the  right  hand  side  of  Eq.  (3.17),  we 
obtain 


T ( t+a)  00  A 

log  -  =  -  (v+i)  log  t  +  £  — 

T(t+a+v+i)  j=l 


(3.18) 


where 


il 


(-1)J 

J(j+1) 


[ B j  +1 ( v+a+x)  -  B j  +1 (a ) ] 


(3.19) 


Thus 


T(t+a)  =  -(v+i) 


T ( t+a+v+i) 


n  +  l  i 

3-1 


(3.20) 


and  can  be  recursively  computed  by  Eq.  (2.7)  as 


]ij  =  1  Ji  1  Ail  Ci0 


(3.21) 


Substituting  Eq.  (3.20)  in  Eq .  (3.17)  and  equating  the  co¬ 

efficient  of  same  powers  of  t,  is  determined  explicitly  as 


1 


(3.22) 


L  Ri-J  =  Qi 


j  =o 


Ro  = 


Now  using  Eq .  (3.17),  (3.15)  in  Eq.  (3.9)  and  noting  that 
the  term  by  term  integration  is  valid  since  a  factorial 
series  is  uniformly  convergent  in  a  half-plane  (Doetsch  (1971)), 
we  have 


f ( u )  *  K (p , n )  K  uSm_1  l  R  2Ti 
1  j  =0  3 

P-1 


d+ic 


d-i  = 


-t  T ( t+a ) 
u  - - 


T ( t+a+v+ j )  (3.23) 


where  =  (2ir)  2  sV  p  P(m+^)+'^>  Now  using  Eq.  (2.4)  in  the 
above  integral,  we  have 


CO 

f (u)  =  K ( p , n )  K  l  R  uSm+a_1(l-u)V+j"1/r(v+j),  0  <  u  <  1. 

■*-  •  —  n  J 


j=o 


(3.24) 


Thus  the  c.d.f  of  u  in  terms  of  incomplete  beta  functions 
lx ( • , ♦ )  is 

»  r ( sm+a ) 

Prob  (u  <_  x)  =  K(p , n)  K..  £  R.  I  (sm+a,v+j)  -  (3.25) 

j=0  3  X  r ( sm+a+v+j ) 


where 


T (a+6) 


a-1 


S-l 


=  ]0  T  (a )  r(B)  u  (1_u)  du  * 


(3.26) 


Further  expansion  of  K(p,n)  = 


f (mp+6p ) 


n  r (m+6-i+l) 
i=l 


and 


T ( sm+a ) 
r ( sm+a+v+j ) 


gives  us  that 


±Z£  .  A* 

log  K  ( p  ,  n )  =  log  [  (  2  tt  )  2  pP(m+l5)"'3  mvj  +  l 

r  =  l  m 


where 


r - 1  1  r+1 

r  r ( r+1 )  i=l  r 


!  I  Br+1(«-1+1)  -  - ] - f 

i  =  l  p  s 


(3.27) 


K(p.n)  =  ( 2 7r )  2  pP(ra+lS)  ^  mV  [1  +  £  -7  ] 


r  =  1  m 


where 


•*:  I  £  *  *  * 

=  r  *  Ai i  ^r-l  ’  ^0  “  1 


(3.28) 


<=°  A , 


rU&tiij)  ■  ^  <«>'<v+3)  *  I  -f  »-2’> 


*  (_1)r  Air 

A.  =  -  [Br+1(a+v+j)  -  Br+1(a)]  =  (3.30) 

sr  r (r+1)  S 


T(sm+a)  ,  ..  N  “  C 

- —  -  (»s)-(',+J)  [1  +  I  -IS] 

r ( sm+a+v+j )  r=l  m 


(3.31) 


1-  c 

*  lr  *  *  ir* 

)r  ■  r  ‘  V  CJ.r-t  ’  s’  10  ' 


(3.32) 


Hence,  Eq.  (3.25)  is  of  the  form 


oo  Q*  oo 

Prob  (u  <  x)  =  (1  +  |  —7)  J  R.  I  (sm+a , v+j ) (ms) 

r =1  mr  1=0  3 


00  v>  , 

1  *  I 

1  r =1  (sm) 


I  (sm+a, v)  +  l  ~ l  Gi 
x  i  =  l  n> 


(3.33) 
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where 


G 

1 


i 

i  I  (sm+a , v+i- j ) 

j  =0  j 


i 


* 


i-Jl 

s 


(3.34) 


The  exact  c.d.f.  can  be  calculated  using  Eq.  (3.33).  When 
n  is  small,  a  suitable  choice  of  5  will  make  m  large  in  order 
to  expedite  the  convergence  of  the  series  in  Eq.  (3.33). 

We  will  now  examine  the  first  few  terms  of  Eq.  (3.33) 
when  n  is  large.  We  know  that 

Co  i  * 

G1  =  I  (sm+a,v)  ( -  +  0,  )  +  R.  I  (sm+a,v  +  l)/s  (3.35) 

X  X  S  XXX 


=  1  (  s  m+  a  ,  v  )  ( - + 

x  i- 


C  o  a  Co i 


Cxi 


+  Q2)  +  R1  Ix ( sm+a , v+1 ) ( — —  + 


+  R2  lx ( sm+a , v+2 ) / s ‘ 


(3.36) 


*  * 

C o  2 

G3  =  Ix(sm+a,v)( 3  +  — 2 + 


Co  3  Q-i  Co  2  Q,Co  i  * 

1  —  +  q3) 

s  s  s' 


C 1 2  Q,Cxi  Q 

+  R^  lx  (  sm+a  ,  v+1 )  (— —  +  - 2 -  +  — ) 

s  s  s 


Cz  i  Q. 


+  R,  Xx ( sm+a , v  +  2 ) (— j—  +  — j)  +  R^  Ix ( sm+a , v+3 ) / s  .  (3.37) 

s  s 

From  Eqs.  (3.27)  and  (3.28),  we  obtain 


*  *  1  „  f 

Q,  =  A  =  -  j  [  l  B,(6-i+l)  - 
1  i=l 


B2(6p) 


(3.38) 


Using  Eqs.  (3.14)  and  (3.16),  we  obtain 


■  A1  =  -sA1 


-  sQ, 


« 


(3.39) 
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Also,  from  Eq.  (3.22),  we  have 


+  Co  1 


(3.40) 


Eqs.  (3.19)  and  (3.21)  give 


Coi  =  Aoi=  [B2(v+a)  -  62(a)]. 

“  2  1 

Using  the  Bernoulli  polynomial  B0(x)  =  x  -  x  +  7-  in 

^  6 

Eqs.  (3.38)  and  (3.41)  and  set 

<$n  *  (2P2+  1 )  / 6p 


We  obtain 


a  =  ( 1 -v ) / 2  =  (3-p  ) /4 
0 

Ai*  -  Ql  “  C01  =  0 


A1  =  Q1  =  R1  =  °‘ 


(3.41) 


(3.42) 

(3.43) 

(3.44) 


Using  Eqs.  (3.14),  (3.16),  (3.19),  (3.21),  (3.22) 
and  (3.44)  we  obtain 

R2  +  C02  =  Q2 


C02  =  A0  2  =  ~6  [  B3  ( v-Ha  0  )  -  B3(a0)] 


(3.45) 

(3.46) 


.2  B,(50p) 


Q2  a2  6  ^  2 


s  r  3  '  “  u  ^ /  _  £ 


i  =  l 


B3(S0-i+l)J.  (3.47) 


By  equating  Eq .  (3.46)  and  Eq.  (3.47)  and  using  B3(x) 


3  3  2  1 

=  x  -  —  x  +  -^  x,  we  obtain. 


,2  =  3P2  (p6-3p4-p2+3) 


4  ( 2p6-6pA+3p2+l) 


(3.48) 


With  the  above  choice  of  s  ,  R2  = 


0. 
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-  2  s  2  =-A2  “  ^2 


where  sfl  is  the  positive  root  of  (3.48).  Hence 


¥  +  Q* 

2  y2 


Now  in  Eq.  (3.22),  we  have 


R3  =  Q3  C03* 

From  Eqs.  (3.14)  and  (3.16),  we  have 

3  *  3  * 

Q3  =  A3  “  “s0  a3  =  -  so  Q3 


’  o  "3 


and  hence 


3  *  3 

r3  *  **s o  (Q3  + 


Using  the  identities  given  above,  it  is  seen  that 


G1  =  G2  =  0 


where 


m  =  n-<5 
0  0 


(3.49) 


(3.50) 


(3.51) 


(3.52) 


(3.53) 


(3.54) 


G_  =  — t  (I  (s  m  +a  ,v  +  3)  -  I  (s  m  +a  ,v))  (3.55) 

J  J  x  0  0  0  xooo 

so 


(3.56) 


Thus  by  so  choosing  6  ,  a  and  s  ,  we  have  the  c.d.f.  of 

00  0 

u  in  asymptotic  form  as 


Prob  (u  <  x)  =  I  (s  m  +a  ,v) 
—  Xooo 


1  R3 

+  — 7  — 7  [I  (sm+a  ,  v+3) 
3  xooo 
mo  s  0 


-  I  (s  m  +a  ,v)]  +  0 (m  )  (3.57) 

^  0  0  0 
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Note  that  when  p  =  2,  then  s  =  1,  a  =  -1/4,  6  =  2/4, 

0  -  0  0 

v  =  3/2.  So,  Eq.  (3.10)  becomes 


r(t+6  )  T ( t+6  -1) 

^  (  t )  =  2  /C  - 2 - o -  (3.58) 

r(2(t+5  )) 

0 

Using  the  duplicating  formula  that 

T(2z)  *  tt-*5  2 2 Z -1  r(z+!s)  (3.59) 


we  have 


4>(  t) 


jj  r(t-(i/4)) 
r(t+(5/4)) 


But,  using  Eqs.  (3.15)  and  (3.17),  we  have 


(3.60) 


♦  (t)  -  ^  21~2<5o  (-1— \  R  ^  (3.61) 

T(t+(5/4))  i=l  1  r(t+(5/4)+i) 

Hence 

•  r(t-(i/4)) 

l  *i  -  -  o- 

i=l  r (t+(5/4)+i) 

For  p  =  2,  Eq.  (3.57)  becomes 


Prob  (u  <.  x)  =  Ix(n-l,3/2)  (3.62) 

which  is  the  exact  distribution  of  w  (see  Nagarsenker  and 
Das  (1975)). 

Krishnaiah,  Lee  and  Chang  (1976)  approximated  the 
distributions  of  certain  powers  of  the  likelihood  ratio  stat¬ 
istics  for  testing  several  hypotheses  in  multivariate  anlaysis 
with  Pearson's  type  I  distribution.  The  accuracy  of  this 
approximation  was  found  to  be  quite  good.  In  the  sequel, 
we  will  refer  to  the  above  approximation  as  KLC  approximation. 
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Tables  1  gives  a  comparison  of  the  accuracy  of  the  approxi¬ 
mations  by  taking  the  first  term  and  the  first  two  terms 
respectively  in  (3.57). 

In  the  table,  the  constant  w  is  defined  by 

c  =  exp(-w/2)  where  P[w£c]  =  a 

and  the  values  of  w  are  taken  from  the  tables  of  Krishnaiah, 
Lee  and  Chang  (1976).  When  the  first  term  in  the  expansion 
(3.57)  is  used,  the  value  of  P(w£c)  is  denoted  by 
whereas  the  corresponding  value  obtained  by  using  the  first 
two  terms  is  denoted  by  a 

TABLE  1 

Significance  Level  Associated  with  the  Asymptotic 
Expression  for  the  Likelihood  Ratio  Test  for  Sphericity 


n 

P  =  3  , 

a=0 . 05 

a=0 . 05 

w 

al 

a2 

w 

“l 

a2 

10 

1.745 

.0499 

.0500 

6.512 

.0476 

.0486 

15 

1.115 

.0500 

.0500 

3.895 

.0495 

.0498 

20 

.820 

.0499 

.0499 

2.791 

.0499 

.0500 

28 

.  576 

.0500 

.0500 

1.924 

.0499 

.0500 
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TABLE  1  (Continued) 


n 

P  =  3 , 

a=0 . 01 

P  =  6 

» 

II 

o 

o 

M 

w 

a  1 

ao 

w 

01 

a  _ 

1 

2 

1 

2 

10 

2.263 

.0100 

.0100 

7.531 

.0092 

.  0095 

15 

1.445 

.0100 

.0100 

4.491 

.0098 

.0099 

20 

1.062 

.0100 

.0100 

3 . 216 

.0099 

.0100 

28 

.  746 

.0100 

.0100 

2.216 

.0100 

.0100 

The  above 

table 

indicates 

that  the 

first  term 

in  the 

asymptotic  expression  alone  gives  a  good  approximation. 

KLC  approximation  is  based  upon  approximating  certain  power 
of  w  with  Pearson's  Type  I  distribution  by  using  empirical 
me mods  whereas  the  approximation  given  in  this  paper  is 


analytic  in  nature 


4.  Asymptotic  Nonnull  Distribution  of  the 
Likelihood  Ratio  Test  Statistic  For  Sphericity 


In  this  section,  we  derive  the  asymptotic  distribution 
of  the  likelihood  ratio  test  statistic  for  sphericity  of 
the  complex  multivariate  normal  distribution  under  the 
following  sequences  of  local  alternative  hypotheses: 

(i)  (I-qI_1)  =  V/m,  (ii)  (I-q-1I)  =  W/m 


where  V  and  W  are  fixed  matrices  as  m-*-00  and  0<q<°°.  The  expressions 
obtained  are  in  terms  of  beta  series.  In  the  analogous 
real  cases,  Khatri  and  Srivastava  (1974)  obtained  asymptotic 
expressions  in  terms  of  chi-square  series.  To  derive  the 
asymptotic  distribution  in  the  complex  case,  we  need  the 
following  lemma  in  the  sequel. 

Lemma  7.  Let  A:  pxp  be  distributed  as  a  complex  Wishart 

til 

CW  (A;n,E).  Then  the  non-null  h  moment  of  w  defined  in 
P 


Eq .  (3.6)  is 


E(wh) 


rp(n) 


n  r  y  C<  T (np+k) 

,4  ^  k !  T (np+k+ph) 

k  =  0  < 


where  M  =  I-qE 

Proof:  The  distribution  of  the  eigenvalues 

of  A  is  (James  (1964)) 


T  (  n+h  ,  < ) 
(4.1) 


K(p,n,E)0F0(-E_1,L) | L [ n_P 


n  U,-*,)2 

i<  j  J 


P 

IT  dJ l  . 


(4.2) 


where 


^(P-D 


r  (n)f  (p) 


K(p,n,I) 
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and 


L  -  d lag  { l  , l  ,  .  .  .  ,  l  ) 

12  p 


Using  Lemma  4,  Eq.  (4.2)  can  be  written  as 

K(p,n,I)  e~tr  L/q  F  (M/q,L) ]L|n_P  H  (l  -l  )2  U  d Z  (4.3) 

i<j  J  i=l  1 


where  M  *  I-qE 

Multiply  Eq.  (4.3)  by  w*1 ,  also  expand  qFq  (M/q  ,L)  as  in  Eq 
(2.9)  and  perform  the  transformation  L  =  UVU '  where  U  is 
unitary  and  V  is  Hermitian  positive  definite  matrix.  The 
Jacobian  of  the  transformation  is 


J ( V ; L , U)  =  n  U ,-l  .)2  h  (U) 
i<  j  2 


and  the  integration  on  U  is 


h2(U)dU  =  irp(p-1)/rp(p) 
U(p) 

(see  Khatri  (1965)). 


We  have 


C«<M> 


1  l 

r  (n)  k  =  0  <  qkCK(I)k:  v4«>0 


exp(-tr  V/q)|v|n+^  p(tr  V) 


X  C|c(V)dV  . 

Applying  Lemma  5  to  the  integral,  we  obtain  Eq .  (4.1).  Note 

that  for  q  =  1 ,  we  get  the  expression  as  in  Pillai  and 
Nagasenker  (1971). 

We  will  first  derive  the  asymptotic  distribution  of 
-1  1  /s 

w  when  (I-qE  )  =  V/m.  Let  u  =  w  5  where  s0  is  the  posi¬ 


tive  square  root  of  the  right  side  of  (3.48)  and  w  was  defined 
by  (3.6)  From  the  non-null  hP^  moment  of  w  given  in  Eq .  (4.1), 
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we  use  the  inverse  Mellin  transform  on  E(u^)  =  E(w^s°), 

to  get  the  density  of  u  as 


°°  C  (M) 

f(u)  =  Kx(p,n,E)  £  l  — -  T(np+k)  *  ^ ( u )  (A. 4) 

k  =  0  k  k  ! 


where 


V  (u) 

and 


2  tt  i 


C+ioo 


c-  ic 


-h-1 


ph/  s 
P  o 


P 

n  r(n+(h/s  )+k  -i+1) 
i- 1 _ 0  1  dh  (4.5) 

r ( np+k+ (ph / s  )) 
o 


K*(p,n,Z)  =  |q  Z  1|n/  if  r(n-i+l)f 

i  =  l 


M  =  I-q  z"1  =  V/m  . 

o 

Following  the  same  argument  as  in  Section  3  we  obtain 


(4.6) 


Prob  (u  <  x)  =  | q  Z 


,  C  (M)  . 

-I  i  n  r  r  <  _  k 


”  Q±(k) 


I  +  I-h--) 

=  0  <  k  l  0  i ~1  m 


* 

«  C  , 


I  *00  I  (s0m0+a  ,v+a)(m0s0)~a[l  +  £  -y1  ]  (4.7) 

a=0  a  X  j-1 


where  we  obtain  the  following  expressions  analogous  to 
(3.22),  (3.16),  (3.14),  (3.28)  and  (3.27): 

X  Ri-J<k)  Ci-j  .  j  ’  qi(k>  ;  Vk>  ■  1 


(4.8) 


Qr(k)  =  7  l  1  A^(k)  Qr-a(k);  Q0(k)  =  1 


(4.9) 


*'  * 


T 
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and 


A  (k) 
r 


("1)r  so  Br+1(60p+k) 


r  (r+1) 


i-1 


Br+l(60+ki“i+1) 3  (4.10) 


Q*(k)  =  r  J  *  AI(k)  Qr-i(k)»  Q0(k)  =  1 


(4.11) 


*  (-nr  ?  Br+1(6ep+k) 

A  (k)  =  - 1—  [  I  B  (6  -i+1)  -  — +^-~ -  ] 

r  r (r+1 )  i-1  r+1  0  pr 


(4.12) 


The  asymptotic  expansion  of  |q  E  3 | n  gives 


q  E_1 ln  =  ( I-V/m  I (m0+60) 


■  exp [- (m  +6  ) (  £ 

n  n 


trV 


o  o 


)] 


r  =  1  rm 


=  exp  (-trV) [ 1  -  —  (5  tr  V  +  J  tr  V2) 

mo  2 


tr  V24tr  V3) 
m  2  J 


-  i  (5  tr  V  +  i  tr  V2)2  } 
2  o  2 


+  0 (m-3 ) 
o 


(4.13) 


and 


C  (M)  m 

K  0 


C  (V/m  ) ra  =  C  (V) 
<oo  K 


(4.14) 


since  the  zonal  polynomial  C^(M)  is  homogeneous  of  degree  k. 
Furthermore  the  use  of  formula  for  the  Bernoulli  polynomial 


B  (x+b)  =  T  r\  B  (x)br 

n  r =0  Vi  n_r 


(4.15) 
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g  ive  s 

Ax(k)  = 


A2(k)  = 


From  eqs . 

Q 


where  A^  , 
(<)  and 
Eq .  (3.44) 


(-l)s  B 0  ( <5  p+k) 

0  r  2  o 


i-1 


B?(6  +k.-i+l)] 
2  o  i 


A  -ys  [~  -  (1  +  i)k  -  a. (<)] 
I  2  o  p  p  1 


(4.16) 


s  B. (6  p+k) 


[- 


3  Q 


- 1 


i=l 


B,(6  +k  -i+1)  ] 
2  o  i 


sj  3  3B  (6  p)k2  3B  (5  p)  ,  . 

S  +  —  +  " — t-0 -  +  ( V"5 -  3  6  -36  -  i)k 

2  6  p2  p2  p2  o  o  2 


-  i  a  2  ( tc  )  -  (36^+  j)  (  k  )  ] 


(4.17) 


(4.11)  and  (4.12) 

Q*(k)  *  A*  (k )  -  A*  +  i  (2  B1(6()p)k+k2)  (4.18) 


* (k)  =  |[A*(k)2  +  2  A* (k) ] 


Q*  +  — (2  B. (5  p)k  +  k2)2 


'2  2 

8p 


'1'  o' 


(4.19) 


- -  (3  B,(5  p)k  +  3  B.  (6  p)k2  +  k3) 

6p 2  2  0  1  0 

A^ ,  Q2  are  as  in  Eq .  (3.14),  (3.27),  (3.28), 
i2(k)  are  defined  in  Eq.  (2.18).  Note  that  in 
,  (3.49) 


A1  =  A1  =  Q2  +  C02  =  Q2  + 


02 


2  2 
s 


(4.20) 


o 


Also  from  Eqs.  (4.8),  (4.9) 


Qx(k)  =  Ax(k) 

Q2(k)  =  j  (A1(k)Q1(k)  +  2A2(k)) 

Rx(k)  =  Q1(k)  -  CQ1  (4.21) 


R2(k)  =  Q2(k)  -  CQ2  -  ( Q 1  ( k )  -  CQ1)  Cll. 


Substitute  these  identities  in  Eq.  (4.7)  and  then  use 

Eqs.  (2.13)  -  (2.17)  for  the  summation  over  k.  By  neglecting 

the  higher  order  terms  of  m  ,  we  obtain  an  asymptotic  ex- 

o 

pression  of  Prob  (u  £  x)  as 
Prob  (u  <_  x) 

*  I  (s  m  -fa  ,v)  +  i—  d-  [I  (s  m  +a  ,\>)-I  (s  m  +a  ,v+l)] 
x  od  o  m  l  x  oo  o  x  oo  o 


3 


l 

i  =  l 


a.  I  (s  m  + a 
i  x  o  o 


v+i-1)  +  0(m 


(4.22) 


where  d,  = 


1  2  12 
(tr  V)  ■  j  trV 


1  1 
a  -  =  6  d.  +  —  +  d  ;  d  =  — 

1  o  1  2  2  2  3p 

,v+l  „„  1  2  tr  V 

a,  =  y—  -  25  -  -  -  - 

2  2  o  0  P  P 


3  1  3 

2  (tr  V)J  -  i  tr  VJ 


>  di  -  di 


«3  ■  <5  +  \  -  iir>  di  +  -7  + 

0 


3p 


(t  >  .  3  +  i  tr  V3  - 


tr  V  tr  V 


1 
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We  now  discuss  the  asymptotic  distribution  of  w  when 
I  -  q  =  W/mQ.  In  this  case,  we  have 


q  I-1  =  (I  -  W/m  )-1  -  I  +  —  (I  -  — )_1 

o  m  m 

o  o 


M  -  !  -  q  r1  -  -  (I  -  ^r1 

mm 

o  o 


(4.23) 


We  now  substitute  Eq.  (4.23)  in  Eq.  (4.4)  and  expand.  It 

may  be  noted  that  when  compared  with  Eq.  (4.6),  -  W(I-WAn  ) 

o 

plays  the  same  role  as  V  in  those  equations  and  hence  the 

non-null  asymptotic  expression  for  the  present  case  can  be 

obtained  by  treating  -  W(I-W/m^)  as  V  in  Eq .  (4.22). 

Further  expanding  -  W(I-W/m  )  3  and  neglecting  the  higher 

o 

order  terms  of  m  gives 

o 

Prob  (u  _<  x) 

=  I  (s  m  +a  ,v)  +  —  d..  [I  (s  m  +a  ,v)  -  I  ( s  m  +a  ,v+l)] 
x  ooo  mlxgoo  X  o  0  0 

0 


-1 


>  J  _  Q 

+  — o  y  a.  I  (s  m  +a  ,v+i-l)  +  0(m  ) 

2  *--  IX  00  0  o 

m  I  =1 


(4.24) 


where  d,  = 


I?  (tr  W)2  -it. 


<  ,  ,1 

<5  d„  +  +  d.  ;  d  -  tr  B  tr  W  -  — , 

o  1  2  2  2  p  3p 


3  2  3 

(tr  W)  J  --  tr  v 


(yr-26 


-  I)d. -d2  +  d_  ;  d_  -  iy(tr  W)3--tr  W  tr  W 
0  pi  1  3  3  p 


‘3  2 

P 


2  3 

+  tr  W 


<1 +  -  i^>  rfi  4 +  ^ 
0 


d,  = 


-  tr  W  tr  W2  -  \  tr  W3  -  — (tr  W)3 
P  3  3p 1 


28 


5.  Applications  in  Inference  on  Multiple 
Time  Series 

In  this  section,  we  discuss  the  applications  of  the 
results  of  Sections  3  and  4  in  drawing  inference  on  the 
spectral  density  matrix  of  the  multiple  time  series. 

Let  x'(t)  =  (X1(t) . xp(t))  (t  -  1.2,.  ...T)  be 

a  real  stationary  Gaussian  multivariate  time  series  with  mean 

I 

vector  0  and  covariance  matrix  R(t)  =  E{X(t)  X  (t+x)}. 

The  spectral  density  matrix  F(u)  at  the  frequency  w  is 
given  by 

00 

F(w)  ~  Yv  £  exp(-ixu>)R(x)  (5.1) 

T  =  -» 

A  well  known  estimate  (e.q.  see  Wahba  (1971))  of  F(w)  is 

A  Si 

F ( to )  =  (f  (oi))  where 

nl 

-  TT^XIT  tj  (5-2> 

1ij(«)  =  ZjL(»)  2  j  (oj) 

1  T 

Z  (u)  = - r  l  X  .  (  t )  exp  ( -itto)  . 

j  ( 2xT) *  t=l  j 

where  Z  denotes  the  conjugate  of  Z.  Then,  it  is  known  (see 
Goodman  (1963)  and  Wahba  (1968))  that  (2n^+l)F(u>)  is  dis¬ 
tributed  approximately  as  central  complex  Wishart  matrix 
with  (2n^+l)  degrees  of  freedom  and  with  E(F(uj))  ■  F(w). 

So,  we  can  test  the  hypothesis  H:  F(w)  ■  A(co)  1^  in  a 
similar  fashion  as  in  the  case  of  the  test  for  sphericity 
of  complex  multivariate  normal  distribution.  So,  analogous 
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to  the  likelihood  ratio  test  for  sphericity,  we  propose 
the  following  procedure  for  testing  H.  We  accept  or 
reject  H  according  as 

T  <  c^  (5.3) 

where 

P[T>Cl|H]  =  (1-ct)  (5. A) 


T  =  F(u)  | _ 

(tr  F(w)/p)P 


(5.5) 


So,  the  approximate  asymptotic  distribution  of  T 


the  null  case  is  obtained  from  (3.57)  by  replacing  n  with 


(2n^+l)  and  is  the  square  root  of  the  right  side  of 
(3.48).  Now,  consider  the  alternative  hypothesis  and 
A2  where  A^  :  ( I-q  ( F  ( w)  )  )  *  V(w)/m  and  A,,  :  (I-q  1F(u)) 

*»  W(ui)/m  where  m  «  2n^  +  1  -  <5  ,  5  is  a  constant,  and 
V(u>)  and  W(to)  are  fixed  matrices  as  m  -*• 


In  the  above  cases,  approximate  asymptotic  expressions 

1/s0 

for  the  distribution  functions  of  T  can  be  obtained  by 

replacing  n  with  (2n^+l)  on  the, right  sides  of  (4.22)  and 
(4.24)  respectively. 

Next  consider  k  frequencies  where  »  2  tt  j  £  /  T 

and  j  ^  =  ( 3.  —  1 )  (2n^+l)  +  (n^+1).  Then,  it  is  known  (e.g.,  see 
Wahba  (1971))  that  F(u^ ),..., F(u>k)  are  distributed  independently 
and  approximately  as  central  complex  Wishart  matrices  with 
(2n^+l)  degrees  of  freedom  and  E(F(tOj))  ■  F(ojj).  We  now 
discuss  procedures  for  testing  and  Hq  simultaneously 


30 


against  and  Ay  where 

H,  :  F(m  )  =  X(w,)I  ,  H. 
i  i  i  p  U 

Ai:  ^  X(ui)Ip,  Aq 

We  accept  or  reject  according  as 


k 

n 

i  =  l 


where 

P[Ti>ci;  i  =  l ,  2 . k|HQ] 

k 

=  n  P[T  >c  I H  ]  =  ( 1 -a )  . 
i  =  l  11 


(5.6) 


_  lFCyl 

1  (tr  F(wi)/p)p 


(5.7) 


The  total  hypothesis  Hy  is  accepted  if  and  only  if  all 
the  individual  hypotheses  are  accepted.  We  can  compute 
approximate  values  of  a  for  given  values  of  c^  using  the 
results  of  Section  3.  The  power  functions  of  the  above 
test  under  certain  local  alternatives  can  be  computed  using 
the  results  of  Section  4.  We  can  also  test  the  hypothesis 
Hy  by  using  Ty  =  as  test  statistic.  Approximate  null 

distribution  of  this  statistic  can  be  obtained  by  following 


the  same  lines  as  in  Section  3. 
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